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QCD phase diagram 


Schematic figure: 


Temperature T 


` Quarkyonic 
Matter 


Baryon Chemical Potential us 
(Fukushima, Hatsuda, Rept. Prog. Phys. 74 (2011) 014001) 


e Finite temperatures --- Lattice simulation works well ! 
e Finite density --- No ab initio calculations. 
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Lefschetz-thimble path integral 


Where is the finite-density QCD? 


Neutron star 

e Cold and dense nuclear matters 

o 2Msun neutron star (2010) 

o Gravitational-wave observations (2017~) 
Heavy-ion collision 


e Low-energy scan of heavy-ion collisions are planned and run in 
many facilities (RHIC, SIS, J-PARC) 


Reliable theoretical approach to equation of state must be developed! 


Z(T, u) = foa Det (P(A) + m) exp —Sym(A). 
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Silver Blaze problem in finite-density QCD 
QCD & |QCD| 


ZQCD = [PA (dety D, ) e 9v. Z|QCD| = [ra Idet^, D, | e ?YM. 


At u = 0, these two are the same! But, 


ng 
QCD 


|QcD| 


m /2 my/3 "H 


At T = 0, the state must be equal to the QCD vac. for |u| S my /3. 
Can we describe this "trivial" behavior using path integral? 
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Motivation 


Path integral with complex weights appear in many important 
physics: 


e Finite-density lattice QCD, 
spin-imbalanced nonrelativistic fermions 


e Gauge theories with topological 0 terms 


o Real-time quantum mechanics 


Oscillatory nature hides many important properties of partition 
functions. 
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Example: Airy integral 
Let's consider a one-dimensional oscillatory integration: 


a= [| ze 1 T as 
p 3 . 


RHS is well defined only if Ima — 0, though Ai(z) is holomorphic. 


Hore T 


exp Os eost) 


AA AAAI 


-8 6 4 2 0 2 4 


Figure: Real parts of integrands for a — 1 (x10) & a — 0.5i 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: (classical eom S'(2,) = 0) 


dx e S@ = n. | dz e 5), 
i 2j, 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


d (=00)| tim. ai) = ca} . det H (8 


No: intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = 24 3). 


[Witten, arXiv:1001.2933, 1009.6032] 
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Lefschetz-thimble integrations 


Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = EL (5 + 2 


Complexify the integration variable: z = x + iy. 


> Re(z) 


4 


Integrand on R, and on A 
(a — 1) 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Oct. 26, 2015 € UConn 9/39 


Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi (5 +a) : 


ng: intersection number of the steepest ascent contour Ke and R. 
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Figure: Lefschetz thimbles J and duals K (a = 1e0-1i, —1) 


E 5 = = = Nae 
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Sign problem & One-site Hubbard model 


Simplest model study to understand the sign problem 
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Sign problem & One-site Hubbard model 


One-site Fermi Hubbard model 


Hubbard model (&: annihilation op., fj = ĉtê): 


= -t o» Creag pe ifia -uY fn; + fry i). 


(1,9) o 


This is too difficult as a toy model. 
Turn off the hopping, t = 0: 


A = Uff, — p(ñ + f). 
Fock state gives the partition function immediately! 


Z=tr [exp -B| — 14 2e8^ + eBQu-U) 
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One-site Fermi Hubbard model 


The number density is given as 


18 Deh + aP a UP) 
B m In Z ui 1 + Jehu + ebQu-U) ` 


In the zero-temperature limit, 


Question 


How can we describe these step functions using the path-integral 
formalism? 
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Sign problem & One-site Hubbard model 


Path integral for one-site model 


The partition function reads 


N p? 
m7 f Teo (32585) 
pr dvd, exp — (>: Viet (vest — ef Mu) : 


with the antiperiodic boundary condition Yy+1 = — 41. 
ọ is an auxiliary field for the number density: 


(A) = (yt, em Get) = — > (pp). 
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Sign problem & One-site Hubbard model 


Integrate out fermions 
Fermion spectrum: 
B N 
Aly, p) = 1 — eG YUN exp 7 Y (gy + w+ U/2). 
k=1 


It couples only to the Matsubara zero mode of y. 


The path integral reduces to the zero-mode integral: 


B : 2 gg. 
Z= do (1 Plip+pu+U/2) Bp? /2U | 
Vo y(l+e Ye 


Integrand has complex phases causing the sign problem. 
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Sign problem & One-site Hubbard model 


Classical action & flow equation 


Effective action: 


S(z)- faan (reos (iie 7). 


It has a charge conjugation symmetry: 


This ensures the manifest reality of observables. 
(YT, Nishimura, Kashiwa, PRD 91, 101701(R) (2015)) 


The flow equation reads 


de £.. 218 exp B (—iz + u + £) 


dt U^ 1+ exp @ (-iz+p+ 9) 
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Sign problem & One-site Hubbard model 


Flows at y/U < —0.5 and u/U > 1/5 


Re(z) 


Figure: Flow at u/U 


2 


| dz e 8), 
Se 


Number density: n, = 0 for y/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 
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Lefschetz-thimble path integral 
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Sign problem & One-site Hubbard model 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at y/U = 0 


-0.25 


-0.5 


-0.75 


Behaviors of the flow become totally different! 


u/U + 1/2. 


~N 


Furthermore, saddle points lie around Im(zm)/U 


This value is far away from n = Im (z)/U=0, 1, or 2. 
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Lefschetz-thimble path integral 
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Sign problem & One-site Hubbard model 


Numerical results 
Results for BU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ja 


r1, and PU J4 


— 1-thimble 
=== 3-thimble 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 


1 UJ 


(cf. Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258) 


Question 
Multiple thimbles can solve the discrepancy. But 
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Semiclassical analysis 


Semiclassical study of the sign problem 
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Semiclassical analysis 


Complex classical solutions 
Let us solve, for —0.5 < u/U < 1.5, 


B. 28 exp B (iz + u+ E) 
iz B 
U 1- exp (iz ++) 


If GU > 1, the solutions are labeled by m € Z: 


U 3U — 
SIE) +T (iem +i Jj + O(T?). 


2 


At these solutions, the classical actions become 
U 12 
Sp x — T (B s 

m ~ 2mm (5 + y f 
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Semiclassical analysis 


Semiclassical partition function 
Using complex classical solutions Zm, let us calculate 


oo 
La t= > p nm. 


Tn—-—oo 


This expression is valid for —1/2 S u/U $ 3/2. 
This is calculable using the elliptic theta function: 


Za 


K 


= 1 dcc 
-S [149 2 H+ \ q-2n?m?/8U 
e + 2, cos 27m (5 + 2 e 


Lo 90 Lad 1 —2n? [BU 
e ^ (r(5+3):. ) 
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Semiclassical analysis 


Number density & Lee-Yang zeros 
Semiclassical study gives the correct transition! 
8 2 (1< u/U < 3/2), 
5; I Za 3 1 (0<p/U <1), 
H 0 CI < u/U < 0). 


Using the infinite-product form of the elliptic theta function, 
=|] {(1— q?5) (1 + 29?" cos(2z) + q^?) . 
(=1 


Zeros of Za. are 


u = +inT(24— 1), p= U +iT(2£=1). 


=> Ist order transition at u = 0, U as T > 0. 
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Semiclassical analysis 


Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za.) = $e | = 6750905 (0, 217), 


o Lee-Yang zeros cannot appear at y/U = 0, 1. 
e One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Ünsal, PRL 115 041601) 
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Semiclassical analysis 


Numerical results 
Let us remind the result for n: 


— 1-thimble 
-=-= 3-thimble 


We need five thimbles to describe the transition at GU = 30. 


Question 
How many thimbles are necessary for practical purpose? | 
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Semiclassical analysis 


Necessary number of complex classical solutions 
Criterion to neglect the Lefschetz thimbles Jm for large |m]: 


m m-—1 
A T 
gm 


cl 


e 


Assuming |m| > 1, we can write this as 


2| exp — Sml <e 
Za ii 


Here e << 1 is a controlling parameter of an error. 
We can solve this inequality with respect to m: 


1 
er Cen) 
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Semiclassical analysis 


Necessary number of complex classical solutions 


If the sign problem is sever, 05 (3, 027/90) ~ yf 288 e-8U/8. 
The criterion gives 


2m 


If the sign problem is mild, 42 ~ 05 (m, 0729/80) pur EE. 


We only need 


Necessary number of thimbles depends on the region which one 
wants to study. 

For most severe cases, it is proportional to BU (— oc). 

(YT, Hidaka, Hayata, 1509.07146) 
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Complex Langevin and Lefschetz thimbles 


Relation between complex Langevin and Lefschetz thimbles 
(to appear soon. In collaboration with Y. Hidaka, and T. Hayata) 
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Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz,(0) OS 


ag 7 gz 00) + V(O). 


0: Stochastic time, 7: Random force satisfying 
(n(8)n(0')), = 20(0 — 0"). 
Ito calculus shows that 


5 (Olen(0)) = A(O" (z,(8))), — (O'(Zn(9))S"(2n(9)))n- 


If the I.h.s becomes zero as 0 > oc, this is nothing but the 
Dyson-Schwinger eq. 
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Complex Langevin and Lefschetz thimbles 


Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(05) 7 7 374. | Mo, 


c o 


ind, € C. To reproduce physics, d, = (Ko, R). 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Complex Langevin and Lefschetz thimbles 


Semiclassical limit 
Let us take h < 1 for computing 


(Ote) = ate [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around zo's imply that 


"ce. st. (0%) S » O. 


RHS is 


[ die SOM"O(z) ~ S Ole). 
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Complex Langevin and Lefschetz thimbles 


Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


C, > 0. And, d, = (Ks, R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, to appear) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points. 
e Those saddle points have different complex phases. 


Consequence 


Naive complex Langevin method cannot explain the Silver Blaze 
phenomenon. 
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Complex Langevin and Lefschetz thimbles 


Proposal for modification 
Assume that 


_ (Ko, R) 2nh So) 


S" (Ze) 


Because of the localization of probability distribution P, it would be 
given as 


P= GPs supp(P,) A supp(P,) = 0. 


E 


Assumption means “CL = phase quenched multi-thimble approx.” : 


(Ko, R) 


27h 
T e SGo)/h 


Sz) O(zo). 


(Olh = 5, 7 


o 
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Complex Langevin and Lefschetz thimbles 


Proposal for modification (conti.) 


If so, defining the phase function 


EME] —iIm S(z0)/h = 
e) (20)/ Xsupp(Po) (2,2), 


we can compute 


new ._ (P(n, 24)O(25))n 
OUT — BG. E. — 


This new one is now consistent within the semiclassical analysis. 
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Possible questions on the proposal 


Our prescription would work well if P is a sum of “well localized” 
distributions and c, can be well understood. 


This prescription has the sign problem because ® is complex, 
but it can be less severe compared with the original one. 
Possible questions are 

e |s the assumption on c, valid? How can we compute it? 


e Can we go beyond the semiclassical analysis? 
If not, what is an error estimate? 


| have no answers.... 


Anyway, let's try! 
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Complex Langevin and Lefschetz thimbles 


Numerical test 1 


S(x) = 24/4 — x? /2 +iax (Semiclassical limit is œ > oo): 
(Preliminary) 


--- exact 
HEH w/o phase 


(Hayata, Hidaka, YT, to appear) 
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Complex Langevin and Lefschetz thimbles 


Numerical test 2 
One-site Fermi Hubbard model: (Preliminary) 


2:5 


--- exact 
20l HEH w/o phase 
tQ w/ phase 


-05 0.0 05 1.0 1.5 


(Hayata, Hidaka, YT, to appear) 


Consequence 


Modified complex Langevin is not perfect yet, but it seems to point a 
correct way. | hope this can attack the Silver Blaze phenomenon. 
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Summary & Conclusion 


Summary and conclusion 
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Summary & Conclusion 


Summary and Conclusion 


e Picard—Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 

e One-site Hubbard model is a nice toy model to play with the 
sign problem. 

e Destructive interference among complex phases of Lefschetz 
thimbles play a pivotal role for the baryon Silver Blaze. 

e Naive complex Langevin method is incorrect for generic cases 
and cannot explain baryon Silver Blaze. 

e Modified complex Langevin may become a multi-thimble 
simulator, which has not been achieved in any other algorithms. 
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